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ABSTRACT 



In recent postings, Foreman & Scott (20111 and Gnedin 
TuUy-Fisher relation (B'l'l^'Rj ot gas rich galaxies as tests o 

These criticisms are rather redundant, as they mostly rehash material I have already discussed. ,Gnedm| 
(2011 1 is concerned with explaining away the problem of the apparently missi ng baryons with ionized 



(2011) criticize my work on the baryoni c 
MDND and ACDM (IMcGaugh 2011a'b" 



gas. This is one hypothetical possibilit y that I have prev i ously discussed ( McGaugh et al _ 
McGaugh fc Wolf|2010| |McGaugh|2011b[ ) . [Foreman fc Scott] ( |2011| claim to find a difference be 



2010 



ween 



the acceleration scale that they rneasure and that 1 measure from the same da ta, but the apparent 
difference results simply from their failure to account for the well known fact ( [Binney &: Tremaine 
19871 that flattened disks rotate faster tha n the equivalent spherical mass distribution (as spelled 
out for this application in McGaugh|[2011b I. They further argue that the the intrinsic scatter in the 
BTFR is merely tiny rather than zero. 'I'his difference is within the uncertainty of the uncertainties, 
so confirms rather than refutes my result. Worse, these papers miss the basic point. Why does MOND 
have any predictions come true? Why does it provide such an effective and economical description 
of galaxy dynamics? They do nothing to address these deeper questions. One can imagine that the 
apparently MONDian behavior of galaxies is merely an approximation to some phenomenology that 
emerges in the context of ACDM, but we have yet to demonstrate how this occurs. Until we confront 
this issue seriously, and confirm the existence of non-baryonic cold dark matter in the laboratory, it 
is reasonable to remain skeptical of ACDM as well as of MOND. 

Subject headings: dark matter — galaxies: kinematics and dynamics — galaxies: dwarf 



1. INTRODUCTION 

In the present cosmological paradigm, ACDM, most 
of the mass-energy of the universe is in invisible compo- 
nents: dark matter and dark energy. Generically speak- 
ing, dark matter could mean any non-luminous gravitat- 
ing mass, baryonic or not. For our cosmological model to 
work, the dark matter must specifically be dynamically 
cold and non-baryonic. 

We have arrived at the ACDM paradigm for very good 
reasons that do not require repetition here. However, to 
date all evidence for dark matter remains astronomical in 
nature, resting fundamentally on the assumption that no 
new gravitational physics occurs when extrapolating the 
force law by eight orders of magnitude from precision so- 
lar system tests to extragalactic systems. The persistent 
presence of mass discrepancies in extragalactic systems 
might therefore also be interpreted as a failure of this 
assumption. 

We come then to a fundamental dichotomy of attitudes 
that different scientists bring to the problem. A common 
one is that ACDM has so many successes that it must be 
correct. Consequently, non-baryonic cold dark matter 
must exist. This is the attitude I brought to the problem. 
However, another valid attitude is the ACDM requires a 
hypothetical new particle generally presumed to exist in 
a hypothetical new dark sector of particle physics. Until 
the existence of these particles is confirmed in the labo- 
ratory, ACDM remains an unconfirmed hypothesis built 
on a good but unproven supposition. Which of these at- 
titudes one brings to the problem inevitably colors one's 
interpretation of any given set of data. 



2. SOME CONTEXT — MY EXPERIENCE WITH THE 
TULLY-FISHER RELATION 

Rotating galaxies obey a re lation betw een their rota- 
tion velocity and luminosity ( Tully fc Fish er 1 977|, or 
more generally, their baryonic mass ( McGaugh et al.| 
2000[ ) of the form 

Mb oc V 



(1) 

with the slope typically found to be in the range 3 < 
a; < 4. Though widely used in distance scale work, the 
physical basis of this relation remains unclear. Starting 
from N ewton's universal law of gravity, [Aaronson et al.| 
( [1979| obtained 

EAf oc t/"* (2) 

where E is the mass surface density. This suffices if, as 
believed at the time, disk g alaxies all shar e essentially 
the same surface brightness ( Freeman|1970 ) so that S ~ 
constant. 

I first encountered this problem in the context of low 
surface brightness (LSB) galaxies, which explicitly vio- 
late the required constancy of S. I had, at the time, de- 
veloped my own (perfectly conventional) theory for the 
formation of LSB galaxies. Motivated by observations 
of the stellar populations of LSB galaxies ( [McGaugh fc 
Bothun|1994 1 that suggested that LSB galaxies were less 
evolved and perhaps had a later start in life than their 
higher surface brightness brethren, I hypothesized that 
their low stellar densities followed from residing in low 
density dark matter halos. Though I arrived at this idea 
independ ently, a similar idea had been previously dis- 
cussed by [Dekel fc Silk[ (|T986| . 

In modern parlance, 1 associated low surface bright- 
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ness galaxies with late forming dark matter halos. This 
hypothesis made two predictions: (1) that LSB galaxies 
should be less strongly clustered than HSB galaxies, and 
(2) that LSB galaxies should be systematically shifted off 
of the TuUy-Fis her relation. Th e first of these hypotheses 
was con firmed (Mo et al .lll994l) and subsequently recon 



firme d (Kosenbaum & Bonians 2004[ Rosenbaum et al. 
2009 1 . I therefore engaged briefly in tne (not uncommon 



fantasy that I understood galaxy formation. 

The data soon falsified prediction 2: LSB galaxies fall 
on the same TuUy-Fisher relation with the same nor- 
malization as higher surface brig htness galaxies (Zwaani 



et al.|199^ Sprayberry et al.|1995i) . Shocked, I struggled^ 
hard to understand fiow this could be, without success 



( [McGaugh fc de Blok||1998| . Every time I thought I 
had succeeded, 1 realized that at some point I had made 
an assumption that guaranteed the desired result. Since 
that time, various papers have been published that assert 
success but appear to suffer from the same shortcoming. 

Though shocked enough at the failure of my own pre- 
diction and the difficulty of recovering a conventional 
model, it was more shocking still to realize that this 
observ ation had been predicted in advance by |Milgroni| 
( 1983 1 . A priori predictions play an important role in the 
philosophy of science because it is often possible to "fix" 
a model once a result is known^. Dark matter models 
are certainly flexible enough to accommodate such fixes. 

Were it permissible to simply ignore Milgrom's predic- 
tion, then I would have patched up my hypothesis and 
claimed success. This would have been convenient, emo- 
tionally fulfilling, positive for my career, but intellectu- 
ally dishonest. I did not get it right a priori. MOND 
did. I relate this story because I value intellectual hon- 
esty very highly. 1 am willing to admit when I am wrong, 
and have done so in the past as I did with my prediction 
for the TuUy-Fisher relation of LSB galaxies. The cur- 
rent situation with the gas rich galaxies is not such an 



occasion. 



3. SPECIFIC CRITICISMS 



Gnedin| ( l2011[ ) and [For eman fc ScottJ pOlTI) make dis 
ict criticisms of iMcUaughi (i20lla|). iGnedinl (|2011|) fo 



tinct criticisms of [McUaughi ( |20lla ). Gnedin ( 2011 1 fo- 
cuses on poten tially unobserved baryons while i^'oremanl 



& Scott (2011 1 focus more on the observed scatter in the 
'l'ully-1 isher relation. I can find little important, factual 
information in either paper that I disagree with. Indeed, 
I find myself in the bizarre position of defending myself 
against points they make as if they were new when I 
have already made many of the same points myself (e.g., 
McGaugh et"al][20l0l [McGau gh 2011b). 

Both seem to misconstrue the meaning of the starting 
equation typically invoked for discussing the BTFR in 
the context of ACDM: 



(3) 



where Ca is an easily derived scaling constant that de- 
pends on the cho i ce of o ver-density A. As I discuss in 
§3.1 of McGaugh (2011b), there does at least appear to 
be a consensus on this relation, which refers to the total 
mass enclosed by A and the circular velocity of a test 

^ Remarkably, Milgrom's prediction was made at a time when 
LSB galaxies were widely thought not to exist because of Freeman's 
Law. 



Mo et al 
of definition m 



particle at the radius enclosing this mass. This is not 
the observed BTFR. 

An obvious and natural route to obtaining a first ap- 
proximation of the BTFR is to define factors fd and fv 
that relate the observed fraction of baryons to the bary- 
onic mass associated with the dark matter halo (M^ = 
fdfbMA) and the observed velocity to Va (^/ = /vVa)- 
While these fudge factors are presumably of order unity, 
they do not have to be identically so. 
This is a standard approach adopted, for example, by 
(1998). The only difference is a trivial one 
they choose A = 200 and subsume 
the universal baryon fraction fb into into their param- 
eter rud = Mb/Af2oo- They assumed rud was constant, 
yielding a BTFR with slope 3. This fit the data for the 
bright galaxies known at the time tolerably well with 
rrid ~ 0.05, somewhat less than = 0.17. This is per- 
fectly admissib le, though it does raise the question that 
Gnedin (2011 ) attempts to address: where are those un- 
observed baryons ? 

The new data for gas rich galax ies implies a much lower 
ru d ^ 0.01 (Gurovic h et al.||2010[ ), as illustrated in Fig. 4 
of McGaugh (2011b|jr In order to match the data, we 
need to make The combination of fd and fv vary with 
velocity as defined in equations 16 and 17 of |McGaugh| 



( 201 lb ) . That is an inference of what is required, nothing 
more. Any model that matches this constraint success- 
fully passes this particular test. 

Indeed, I constructed such a mod el myself over a 
decade ago. Starting from the model of Mo et al. ( 1998 ), 
I noted that the BTFR could be fit by inaking rud a lin- 
ear function of circular velocity. I never published this 
model because it is obviously wrong. Though one can in- 
deed fit the BTFR by fine-tuning nid oc Vf, this destroys 
other virtues of the model. For example, the disks size- 
circular velocity relation, which is well explained with 
constant m^, now fails (Fig. [T]) because the dependence 
of disk scale length Rd on rrid and circular velocity in 
their equation 12 cancel out. 

I have spent a great deal of time over the past 17 years 
attempting to construct acceptable ACDM models. I 
have not succeeded. There are no lack of those who have 
claimed success, but these models are rarely consistent 
with one another. At most one of these models can be 
correct, so it seems more likely that the claims o f success 
are overstated (e.g., Trujillo-Gomez et al.||2010 ) 



3.1. Gnedin 



Gnedin (2011 1 does not claim to present a model that 
matches the data. He merely argues that plausible as- 
trophysical processes will reduce /^^ below unity. I agree 
entirely that this is possible (Fig. 12]). If so, it is incum- 
bent upon us to find and identify those baryons, not 
merely accept that they are there is some unobserved 
form^ (in this case, ionized gas). This was explicitly dis 



cussed in [M cGaugh et al.|(|2010D. Indeed, iMcGaugh fc 

■ im 



Wolf ( |201D ) have already published (their t ig. 4) a ver- 
sion of Gnedin's Fig. 2 that shows that re-ionization alone 

^ The ACDM missing baryon problem is not subtle. For every 
detected baryon in a ^ 20 kms"'^ galaxy, ~ 50 remain undetected 
iMcGaugh et al. 2010). Much of the community seems to accept 
tnis without blmkmg. Yet when MOND suffers a facto r of ~ 2 
missing baryon pro blem — as it does in clusters of galaxies l|Sanders| 
|fc McGaugh|2002[ | — we are quick to believe it is falsifiecf! 
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is inadequate to explain the variation of fd- This is one 
example of what I mean to find it bizarre to address as 
criticisms points I have already published myself. 

3.2. Foreman & Scott 



The criticism of Foreman & Scott (2011) is consider- 
ably less coUegial m tone. I'hey appear to be suffering 
from an ancient emotional impulse: when you don't like 
the message, kill the messenger. Anyone who has at- 
tempted to discuss the evidence for evolution or global 
warming in a public setting has no doubt experienced 
this impulse first hand. 

A classic approach that is employed when one wants 
to belittle an argument is to focus on a tiny piece of 
it, show it is wrong, and thus cast doubt on the en- 



tire work. This is the approach Foreman & Scott (2011 1 
adopt. They analyze the the same data as myself'^ and 
find that the peak Vf/GMb = 1.59 As'^ is different 

1.3 As 



Mc- 



from the gp = 1.3 As ^ that I find (see Table 2 of 
Gaugh[2011bp. This is not some horrible mistake as they 
wbUlcTapparcntly like to believe. As discussed in §3.2.2 
of jMcGaugh (2011b ), a paper they cite but seem not to 
have read, the acceleration scale estimated in this man- 
ner is identical to the MOND acceleration scale only in 
the limit of infinite distance from an isolated point mass. 
Since such a condition in never obtained in reality, there 
is a correction factor of order unity that follows sim- 
ply form the fact that a flattened mass distribution like 
a disk galaxy spins faster than the equivalent spherical 
mass distribution. This is a basic result that should be 
familiar to any student of galactic dynamics ( Binney & 



Tremaine||1987 ) 

tween the acce 



This factor ( 0.8) is the difference be- 



eration scale of Foreman & Scott (20111 



a nd myself. 



Foreman & Scott 



J 2011) also argue that the scatter 
ightly larger than I find. While I 



certainly consider the problem of the missing baryons to 
be a puzzle, it is at least approachable. The greater prob- 
lem for ACDM, to my mind, is the small intrinsic scatter. 
Once we allow for undetected baryons {fd < 1), as we 
must in order to reconcile the ACDM halo mass-rotation 
velocity relation with the observed BTFR, then any part 
of the TuUy-Fisher plane below the line Mb = fbCAV^ 
is accessible to ACDM models (Fig. [2]) . It is easy to con- 
ceive of models where galaxies of a given circular velocity 
have fd ranging from near unity to near zero. Indeed, it 
is not uncommon to see such a range in numerical sim- 
ulations. The scatter of the resultant BTFR would be 
measured in dex rather than tenths of a dex. So the real 
question is why the scatter is so small. How do galaxies 
"know" to accrete just the right fraction of cold baryons 
for their rotation velocity? 

Though I am hardly the first to find negligible scatter 
in the Tully-Fisher relation (e.g., Vcrheijen 2001; Ander- 



sen et al 



2001 1 Trachternach et al.||2009| ), ilbreman fcj 



T) make the case that there is some scatter m 



the data 1 used: 0.28 dex rather than the 0.24 dex ex- 
pected purely from the median observational error. Tak- 
ing this at face value, and the difference in quadrature, 
leaves an intrinsic scatter of 0.15 dex — precisely the 



value I give in McGaugh (201 lb I. They also analyze a 



larger sample that has a somewhat larger scatter — 0.33 
dex. If the typical errors are no larger, this implies an in- 
trinsic scatter of 0.22 dex: still a remarkably small value 
given the potential variation in fd alone (neve r mind the 
many other potential sources of scatter, e.g., |Eisenstein| 



fc Loeb._1996) . 
1 do not know the contents of the larger sample [Fore- 



man & Scott ( 201l| ), so I cannot repeat their analysis. 
However, if it is indeed composed of gas rich galaxies, I 
would surmise that a goodly number of the additional 
galaxi es are likely to come from the data of Begum et al.| 
(20081. These galaxies are wonderful for being some of 
the most extreme (in terms of low velocity) yet observed. 
They are also very difficult observational targets. One of 
the harder parameters to adequately constrain is the in- 
clination, which can easily lead to systematic uncertain- 
ties that invalidate the usual assumptions of data anal- 
ysis. A s discussed at considerable length in McGaugh] 



(2011b), I have analyzed only galaxies for which 1 am 



www.astro.umd.edu/~ssm/data/gasrichdatatable.txt 



reasonably confident of the inclination estimate, as indi- 
cated by broad agreement between optical and HI data. 
If one includes all of these data, the scatter certainly goes 
up. The best-fit s lope also becomes somewhat shallower 
( McGaugh|201 lb ) , as these data skew to low velocity for 
a given mass. I'his skew is the expected signature of sys- 
tematic inclination errors, as even a face-on galaxy may 
not appear perfectly circular if it is affected by an oval 
distortion or bar. It is hard to know for certain whether 
this is what is going on, but it is hardly surprising that 
the scatter goes up incrementally as progressively less 
reliable data are included. 

Dickering about whether the intrinsic scatter is zero or 
~ 0.2 dex is of little consequence when variation in fd 
could easily lead to an order of magnitude larger scatter. 
As it is, we have no clear prediction for the BTFR in 
ACDM, let alone the scatter therein. One can certainly 
build models with little scatter, but one must be care- 
ful (as I have repeatedly found) that we do not make 
an assumption that builds-in the desired result. Histori- 
cally, this has routinely been done by simply assuming a 
constant value for fd- 

One common approach to explain fd<lis to invoke 
feedback from supernovae. This is not well understood 
and can have many different implementations. The ques- 
tion then becomes how well these implementations rep- 
resent reality. I would expect a messy process like su- 
pernova feedback to result in considerable scatter — a 
very small gas rich galaxy might experience zero, one, 
or numerous supernovae. Such differences should propa- 
gate into scatter in the BTFR, where the limited scatter 
budget can be consumed very rapidly. 

^Foreman fc Scot t ( 2011) also claim to find higher 
than 1 do. 1 do not see an error in my code. I do notice 
that if I fail to account for the slope of the relation in the 
computation of the weights on the bivariate errors, then 
I get the same number for that they get for the same 
sample. 

For the sake of argument, let us suppose that the higher 
xt ~ 1.5 is correct. Now suppose that it is ACDM rather 
than MOND that predicts a relation with no scatter. Do 
we claim that this falsifies ACDM? Or do we report it 
as a success for the paradigm, with maybe a hint that 
there is some small contribution to the error budget that 
hasn't been entirely accounted for? 
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I have gone through the statistics of these data, and 



the var ious subsamples, in excruciating detail (McGaugh 
2011b). Doing so again would not prove anything, so 1 
make use of an entirely indepen dent program written by 
Ben Weiner ( Weiner et al. 2006 1 . Treating all parameters 
as free gives a best fit slope of 3.82 ± 0.22 with ^ = 
33.5 for 47 galaxies. This slope is consistent with what 
I find myself for the same fit parameters for the same 
data, but is somewhat lower than I found. Fixing the 
normalization and slope to the MOND values yi elds = 
58.5 — sli ghtly more than I found but less than [Foreman 



fc Scott I (120111) find. While it is annoying that there 
should be even a slight disagreement in such a simple 
calculation, the difference is within the uncertainty of 
the uncertainties. 

It is a well known fact that in astronomy, some data 
are more r eliable than others . Rat her than expanding the 
sample as Foreman & Scott ( 2011| ) do, one needs, if any- 
thing, to restrict i t further Doing so essentially returns 



us to the results of Stark et al. ( 2009 1 , who held the data 
to the highest standards ot quahty control. Their empir- 
ical resuh is consistent with MOND'' {x = 3.94 ± 0.07). 
In trying to extend the sample , I have include d some, 
but not all, of the data from Begum et al. (20081. These 
are difficult observational targets th at are inevitably les s 



reliable than the data discussed by Stark et al. (|2009 ) 



Indee d, only three of the twenty-nine galaxie s ot |Begum 



et al. (12008 ) meet the quality requirements of Stark et al 



IgOgf nS discussed in [McGaugh| ( |20 1 lb| , I reject about 



half the sample as having unreliable inclinations. An 
example of just how difficult it can be to measure the 
inclination of these objects is given by the case of Holm- 
berg II (Fig. [3| , which has rather better data than the 
majority of comparably low mass galaxies. How ever , in- 



clination es timates for it vary from 31° to 45 ° (|de Blok 
eral]|2008| to 55 
Carignanl 12002 



' or even 84° at large radii (Bureau &: 
002[ ). I do not make use of data with such 
large disparities in their inclination estimat es. In order 
to expand th e size of the sample as they do, [Foreman fc:| 



vo ^ ^ 

Scottl (2011 ) must invariably do so. It therefore comes 



as no surprise that they find a larger scatter and higher 
, as these are the expected results of incorporating less 
accurate data. 

Returning to the issue of the s catter, one of the brilliant 
aspects of the program of Weiner et al. (2006) is that 
it optionally allows one to ht for the intrinsic scatter. 
Using this option returns a best fit intrinsic scatter of 
zero. Zero is the preferred amount of intrinsic scatter 
the program returns irrespective of whether I restrict the 
slope and normalization to the MOND values or leave 
them free. Weiner's program provides no independent 
substantiation of the claim of finite intrinsic scatter made 



b y [Foreman fc Scott] ( |20Tll ) . 

1 value intellectual honesty very highly. If I had found 
any serious error in my work, I would be eager to re- 
solve the discrepancy and say so. Perhaps I've missed 
something, but I do not see what. Stepping back from 
the details. Foreman & Scott (2011) find essentially the 



This is a good example of wanting — or in this case, not 
wanting — a particular result. I did not want to find a slope of 4 
when we did this experiment because I knew it would eventually 
lead to exactly this sort of dispute. The data didn't give us liberty 
to say otherwise, however. 



same result as I do: an acceleration scale consistent with 
MOND with very little intrinsic scatter in the BTFR. 
They just don't like it. 

4. THE REAL POINT 

I am happy to consider ACDM models that ma ke a 
legitimate effort to expla in the data. Gnedin (2011) and 



Foreman & Scott (2011) fall short of this standard. In- 
deed, they do no even attempt to meet it. 

The point is that MOND provides an economical de- 
scription of galaxy scale kinematics. The simple formula 
proposed by Milgrom is the effective force law in rotating 
galaxies. This is what we need to explain, in ACDM or 
an other theory. 

It is not easy to understand how this simple formula 
follows from the complex physics of galaxy formation in 
ACDM. If it were, there would already be a well estab- 
lished model for it. Instead, the literature is littered with 
mutually inconsistent models that persistently fail to de- 
scribe galaxy data as well as MOND. 

The BTFR connects the observed baryonic mass to 
the characteristic circular velocity. This relation between 
global quantities is just one manifestation of MONDian 
phenomenology. There is also a point by point mapping 
between the observed distribution of baryons and the ob- 
served rotation curve, including all the bumps and wig- 
gles in both. This is a reiteration of the well established 
observational fact that MOND fits the rotation curves of 



galaxies (Sanders & McGaugh 2002 



One can write this as a scaling relation v from the ro- 
tation curve that Newtonian gravity predicts for the ob- 
served baryonic components Vb(r) to the actual observed 
rotation curve Vc(r): 



(4) 



As I've written it he re (see, e.g., Milgrom fc Sanders] 
[2008'; McGaugh '2008), gN = V^/r is the Newtonian ac- 
celeration produced by the observed stars and gas as de- 
termined by numerical solution of the Poisson equation. 
There is no approximation (e.g., exponential disks); the 
actual observed surface density is used, including bulge, 
disk, and gas, with all their bumps and wiggles. For rea- 
sons not presently understood, the simple, smooth func- 
tion V maps from the observed baryon distribution to the 
observed rotation curve. Only the observed baryons are 
required as input, with no reference to the dynamically 
dominant dark matter. 

Stated this way, equation [4] is merely a restatement of 
MOND. However, the scaling relation still holds even if 
MOND is incorrect, as it is known to fit the majority 
of rotation curves. There is only one tiny bit of leeway. 
Computing the baryonic rotation curve 



V^{r)^Vhr) + V^{r) 



(5) 



requires knowledge of the mass-to-light ratio of the stars. 
(Note that this does not apply to the gas, but one can in 
principle have separate mass-to-light ratios for the bulge 
and disk.) This stellar mass-to-light ratio T* is the one 
fit parameter of MOND rotation curve fits. Since we 
are now considering the case where such fits are just a 
scaling relation, not a fundamental law of nature, there is 
no need for this to be the correct mass-to-light ratio. We 
can generalize the relation to account for this by defining 
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the ratio 



Q = 



yTRUE 
yMOND 



(6) 



which is the ratio of the actual stellar mass-to-light ratio 
to that required to obtain a MOND fit. 
We can now write 



(7) 



At least I am trying to address the problem. My critics 
are simply ducking it. 

It is easy to believe that MOND is incorrect as a funda- 
mental theory. I have reviewed many dozens and perhaps 
hundreds of claims to falsify it. Most are unconvincing 
or obviously wrong. Two hold water: th e resid ual mass 
discrepancy in clusters of galaxies (e.g., Sand ers 2003 

- the bullet cluster ( [Ciowe et a 



2007 



Angus et al. 2008 ) 



I 



This mass-discrepancy-acceler ation relation is purely 

It is 



empirical and generally valid (McGaugh 
equivalent to MOND in the case ot y 



2om. 

but it holds 



2004 1 is merely a dramatic example of a more general rul e 
— and the third peak of the CMB ( |Spergel et al.||2007[ ). 
It might be possible to explain these m MUINU, just as 
it may be possible to explain MONDian phenomenology 
in ACDM, but ACDM is the more natural explanation 
for these particular observations. 

Do clusters and the CMB suffice to prove the exis- 
tence of non-baryonic cold dark matter? Certainly not 
clusters — there is nothing about the missing mass in 
clusters that specifically requires the unseen mass to be 
non-baryonic^. The third peak of the CMB is more per- 
suasive, as it implies the oscillation of a mass component 
not coupled to the baryons. That does say cold dark 
matter to me, but I cannot claim to rule out the possi- 
bility of similar behavior from some additional scalar or 
tensor field that might be introduced to combine MOND 
with General Relativity ( Bekenstein||2004 Skordis et al. 



20061 



Definitive proof of the existence of non-baryonic cold 
dark matter requires laborat ory detections. Rece ntly, 
the CRESST-II collaboration fA.ngl oher et al ||2011| ) has 
claimed a detection of WIMPs. Unfortunately, there is 
as yet no corroboratio n: their claimed det ection violates 
limits set by CDMS-II jAhmed et al.|2010|), EDELWEISS 



(Armengaud et al. 



rrr 

UKMA (e.g., [Savage et al.||2009i ). Clearly there remains 



2011), and XENOMIOO (Aprile et al. 



even when MOND does not. This encapsulates the fact 
that it is usually possible to obtain a reasonable fit to 
rotation curves with MOND. The MOND mass-to-hght 
ratios are fairly r easonable in terms of stellar populations 
( McGaugh]|2004 ) , so presumably Q is not very different 
from unity, though it could be. Note that this empiri- 
cal formulation need not apply to non-rotating systems; 
perhaps it is a scaling relation specific to disk galaxies. 

While we need not explain MOND in ACDM, we do 
need to explain the data. The data for disk galaxies are 
encapsulated by equation [7] A satisfactory model would 
show how this equation emerges from the complicated 
process of galaxy formation. It need not even hold all the 
time — there could be some intrinsic scatter around this 
relation. But it does need to hold very nearly true most 
of the time if we are to achieve an acceptable description 
of observational reality. 

We are far from realizing this ideal. In order to apply 
this test to a ACDM model, said model must make a 

detailed prediction for the distribution of baryons from 

which Vg{r) and V^,{r) are calculated. Yet the distri- while b eing difterent from the claimed detection of' 

bution of baryons is notoriously difficult to predict in 
numerical simulations. Until we have models that make 
such a prediction and actually pass this test, it seems 
rather sanguine to assume that it will all magically work 
out. 

Indeed, what disturbs me is the degree of fine-tuning 
the observed phenomenology implies. MOND is the ef- 
fective force law in disk galaxies. The inverse-square law 
is the force-law in the solar system. If we concoct some 
theory of gravity that is different in the solar system, say 
an inverse-cube law, then we would be obliged to say that 
the effective force law only happens to look like an inverse 
square law because there is dark matter arranged just so 
that it always appears to be the case. We would reject 
such a hypothesis out of hand as unacceptably fine-tuned, 
yet the situation is analogous in galaxies. We know the 
effective force law, and must perforce arrange the dark 
matter just so that we can always predict the rotation 
curve from the detailed distribution of baryons alone. 

I have struggled hard to understand how MOND phe- 
nomenology can arise in ACDM. So far, I have failed^. 

^ Few others have even tried. I am only aware of three efforts in 
the literature, each less convincing than the previous. Motivated 
by a p reliminary version of the mass discrepancy— acce l eratio n re- 
lation jMcGaugh 1999), v an den Bosch fc Dalcantonl ( |2000| l ran 
simulations until they tound a teedbaclc prescription tnat leaves in- 
tact only those galaxies that do not deviate too much from it. If 
they are correct, then all other models should use the same feed- 
back prescription. They assumed purely exponential disks, so the 
ability of the MOND formula to fit bulg es and other bumps and 
wiggles remains undemonstrated. |Kaplinghat fc Turner] l |2002| as- 



a lot to sort ouf 

If WIMPs are detected in the laboratory, then an essen- 
tial tenant of ACDM is confirmed. In this case, we still 
need to explain the phenomenology en capsulated in th e 
mass discrepancy-acceleration relation ( McGaugh|2004 ) . 
I desperately hope that we detect WIMPs soon so that 
we'll know and can get on with it. 

Of course, it is dangerous to want a particular result 
in science. What if we do not detect WIMPs? At what 
point do they become untenable as dark matter candi- 
dates? If we are forced to abandon WIMPs, what then? 
Axions? If those don't pan out, do we just invent some 
other, even harder to detect particle? If that fails, why 
not do it again? And again and again, ad infinituml 
While it is possible to exclude specific dark matter can- 
didates, the concept of dark matter generally is not fal- 
sifiable. If it is not falsifiable, does it qualify as science? 



sert that the acceleration scale of MOND follows naturally from 
ACDM, arguing that the acceleration caused by galaxy scale dark 
matter halos is a]j]\j ~ cHq (their equation 4). This does nothing 
to demonstrate the much stricter relation of equation [7] It also 
misses the essential point that the MOND sc ale is connecte d to 
a critical baryonic surface density (ao = GS^ : [Milgrom [198911 for 
which there is at present no explanation in ACJJM. |Dunkel| | |2004[ l 
appears simply to present a tautology. 

® Apparently it needs to be collision-less in the bullet cluster; 
baryonic MACHOs would satisfy that. The collision velocity of 
the bullet cluster (Angus & McGaugh 2008) and t he morphol- 
ogy o ther bullet-like systems (e.g., I'andora's cluster: ,Merten et al.| 
|2011[ ) complicate the interpretation for any paradigm. 
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Figure 1. The disk scale length-rotation velocity relation (analogous to Fig. 4 of |Mo et al.|199"8l l. Models of they type |Mo et al.| ( |1998^ 
discuss do a good job of explaining the data. Both the mean trend and the scatter are well matched by a constant baryon fraction rrifj (sohd 
orange line) with the range of spin parameters A expected from numerical simulations (dashed orange lines). However, in order to fit the 
BTFR, one needs to drop the assumption of constant and scale it as oc Vf. This fits the BTFR, but destroys the consistency with 
the scale length-rotation velocity relation (horizontal line. The normalization of this line is adjustable but not the slope.) This illustrates 
the general principle that one must be careful that tuning done to fit one observation does not upset consistency with another. 
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Figure 2. The BTFR with all data discussed by |McGaugh| l [2011b[ |, making no distinction between gas and star dominated galaxies. 
Galaxies identified there as potentially having large systematic mclmation errors are shown as smaller symbols. Inclusion of these objects 
will increase the scatter and x^i and may also skew the slope. Also shown is the ACDM mass-rotation velocity relation (straight line). This 
shows all of the baryons enclosed within the overdensity A = 200. ACDM does not require galaxies to fall on this line, though that would 
have been a very satisfactory result. If not all baryons are detected (/^ < 1), galaxies may fall anywhere in the shaded region. That they 
fall along such a narrow track when so much parameter space is available strikes me as a fine tuning problem. This is only compounded by 
the fact that the data follow the particular BTFR, that is unique to MOND. The dark gray region represents the allowed regime according 
to Gnedin (2011) — the curved line is his equation 1 with q = 1. Reionization may well reduce the cold baryons available to be detected 
in low mass halos, but by itself it does nothing to explain the observed BTFR. Note that if Gnedinj (j2011jl is correct that this is an upper 
limit to the baryon content of galaxies in ACDM, the two galaxies left of the curved line (which 1 independently identified as having very 
uncertain inclinations) are impossible in A CDM. Indeed, the line moves right and through the data as a increases, so a > 1 can be excluded 
(see also |Tikhonov fc Klypin|2009||McGaugh fc Wolf|2010| l. 



Figure 3 . The d warf galaxy Holmberg II in stars as traced by the Kg image from 2MASS | |Jarrettj ;t al. "200 31 left) a nd atomic gas from 
THINGS ( [Walter et al. 2008 right). The inclination of this galaxy is listed variously as 31° or 45" Ide Blok et 3. |2008| l or more ( see text). 
I have refrained from utilizing data for cases where the inclination is this unsettled (McGaugh 2011bJ; the larger scatter and |Foreman| 
|fc Sco ttl l|2011| find probably stems from including such cases. In this particular galaxy, larger inclinations pose a problem for MUINU 
llbanchcz- y alccdo & Hidalgo-Gamcz 2011). Consistency with the BTFR requires i 25° (green ellipse). That is, to believe that this galaxy 
poses a challenge to iViUlN U requires you to believe that the inclination of this galaxy is significantly different from that illustrated by the 
ellipse, which is not a fit, just the projection of a circle inclined by 25°. Note that masses in MOND scales as [Vobs/ sin(i)]'*, so even a tiny 
error in inclination can lead to a large error in mass for low inclination galaxies. 



